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基于 修正 偶 应 力 理论 和 微分 变换 法 
的 输 流 微 管 自由 振动 
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摘 ， 要 :考虑 黏 性 流体 在 微 管 道内 作 层 流 运动 ,给 出 了 黏 性 流体 在 微 圆 管 中 的 速度 分 布 方程 。 利 用 
修正 偶 应 力 理论 和 Euler 梁 模 型 建立 细 长 微 管 模型 ,根据 虚 功 原理 推导 输 流 微 管 流 - ERO ARRA 
程 , 应 用 微分 变换 法 计算 微 管道 系统 的 固有 频率 。 通 过 与 有 限 差分 法 求解 结果 对 比 ,证 明 微 分 变换 
法 具有 较 高 的 精度 。 随 后 ,研究 了 流体 夭 性 、 微 管材 料 内 豪 特 征 尺 寸 和 预 应 力 对 国有 频率 的 影响 。 
最 后 ,分 析 了 流体 临界 流速 与 预 应 力 的 关系 。 数 值 结 果 表 明 : 在 流体 平均 速度 相同 的 条 件 下 ,考虑 
流体 黏 性 时 微 管 各 阶 固 有 频率 偏 低 , 并 且 平 均 速度 越 大 ,这 一 趋势 越 明 显 。 
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Free vibration of fluid-conveying micro-tubules using the modified 
couple stress theory and differential transformation method 
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Abstract : Considering a laminar flow of viscous fluid in micro-tubules , the velocity distribution equation of 
viscous fluid in a circular micro-tubule is provided. The modified couple stress theory and Euler beam 
model are used to develop the slender micro-tubule model , and the liquid-solid coupled vibration equation 
of the flow-conveying micro-tubules is derived based on the principle of virtual work. The natural frequen- 
cies of the flow-conveying micro-tubules are calculated using the differential transformation method. It is 
indicated that the differential transformation method is of high accuracy compared with the results obtained 
by the finite difference method. Subsequently , the effects of liquid viscosity , material length scale parameter 
and pre-stress on natural frequencies are examined. Finally , the relationship between the critical speed of 
fluids and the pre-stress is analyzed. From the numerical results ,it can be concluded that under the condi- 


tion of the same average velocity of fluids, the natural frequencies of micro-tubules become lower when 


收 稿 日 期 :2022-04-19 修 回 日 期 :2022-09-07 
基金 项 目 :国家 自然 科学 基金 资助 项 目 (No. 11972240) ;广西 密码 学 与 信息 安全 重点 实验 室 研 究 课题 ( No. GCIS201905 ) 
通信 作者 ; 李 成 ,教授 。E-mail: licheng@ czust. edu. cn 
引用 格式 : 朱 成 秀 , 随 岁 寒 , 李 成 . 基于 修正 偶 应 力 理论 和 微分 变换 法 的 输 流 微 管 自由 振动 [中 .应 用 力学 学 报 ,2023 ,40(2) :443-449. 
ZHU Chengxiu, SUI Suihan, LI Cheng. Free vibration of fluid — conveying micro — tubules using the modified couple stress theory and differ- 
ential transformation method[ J ]. Chinese journal of applied mechanics ,2023 ,40(2 ) :443-449. 


nds 应 用 力学 学 报 


NO 
LninaxIV 五 1 


considering the fluid viscosity ,and the higher the average velocity is ,the more obvious the downward trend 


of natural frequencies is. 
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近年 来 ,有 大 量 文献 研究 输 流 管道 振动 与 稳定 
性 问题 5?1 ,其 中 也 包括 微 纳 米 输 流 管道 耦合 振动 
的 研究 。 文 献 [1] 研究 宏观 悬臂 输 流 管 的 动力 
学 行为 ,利用 Galerkin 方法 求解 运动 方程 ,分 析 了 脉 
动 内 流 激 励 与 基础 激励 联合 作用 对 悬臂 输 流 管 的 影 
响 。 本 研究 的 输 流 微 管 可 视 为 人 体 血管 动力 学 的 基 
本 模型 。 此 外 , 随 着 微 加 工 技术 的 发 展 ,在 热 交 换 器 
领域 也 有 着 广泛 的 应 用 ,由 于 输 流 管 的 微型 化 ,导致 
微 管 流动 基本 是 层 流 。 研 究 输 流 微 管 的 力学 行为 对 
币 系 统 的 设计 有 着 重要 意义 。 文 献 [8] 以 输送 无 攻 
不 责 压 缩 流体 的 固 支 单 层 碳 纳 米 管 为 研究 对 象 , 考 
虞 温度 场 和 轴 向 磁场 作用 ,应 用 非 局 部 理论 分 析 模 
型 稳定 性 。 文 献 [12] 考 虑 非 局 部 理论 和 表面 效应 ， 
建立 了 一 种 新 型 纳米 输 流 管 自由 振动 解析 模型 ,分 
析 纳 米 管内 外 表面 层 的 表面 弹性 和 残余 表面 张力 对 
自由 振动 的 影响 。 文 献 [ 13] 建立 了 悬臂 纳米 输 流 
管 术 型 ,增加 端 部 的 从 动力 以 及 切 向 分 布 载荷 作 用 ， 
老 虞 非 局 部 弹性 理论 ,分 析 输 流 管 临 界 流速 的 影响 
四 潮 。 不 同 于 以 上 分 析 纳 米 管 等 纳米 尺寸 量 级 的 输 
流 管 , 本 研究 对 象 微 管 属于 微米 尺寸 量 级 , 主要 基于 
血管 以 及 微型 热 交 换 器 中 输 流 微 管 的 实际 尺寸 。 
(已 有 输 流 管道 振动 文献 多 以 固体 力学 研究 为 
了 通过 假设 流体 无 黏 性 , 并 在 管道 中 匀速 流动 ,以 
简 腰 流体 的 影响 。 实 际 上 , 若 考虑 流体 黏 性 ， 
流 笨 在 圆 管 中 层 流 的 速度 分 布 是 以 圆 管 的 轴线 为 中 
心 线 的 二 次 抛物 面 , 这 和 无 黏 假设 有 很 大 区 别 。 为 
更 加 贴 合 工程 实际 ,有 必要 在 流 固 耦合 振动 的 研究 
中 考虑 实际 流体 的 黏 性 因素 。 针 对 黏 弹性 流体 , 文 
献 [15] 建 立 输送 黏 弹性 流体 的 灾 曲 纳米 管 的 耦合 
模型 ,研究 非 局 部 参数 .曲率 以 及 流体 竺 性 对 模型 振 
动 的 影响 。 文 献 [16 ] 研究 在 热 环境 下 承受 轴 向 谐 
波 载荷 并 输送 黏 性 流体 的 纳米 管 的 动态 稳定 性 , 导 
出 了 耦合 运动 方程 并 分 析 了 各 参数 的 影响 。 这 类 文 
献 多 基于 非 局 部 理论 推导 运动 方程 ,其 中 非 局 部 
参数 在 各 种 不 同 模型 和 工 况 下 往往 难以 定量 给 出 。 
本 研究 将 输 流 微 管 简 化 为 Euler 梁 模 型 ,基于 修正 偶 
应 力 理论 ,研究 输 流 微 管 的 振动 特性 。 利 用 虚 功 原理 
导出 夭 性 流体 作用 下 的 输 流 微 管 系统 控制 方程 。 
在 求解 方法 上 ,以 往 对 输 流 管道 振动 问题 主要 


采取 有 限 元 法 \、 微 分 求 积 法 “等 ,而 微分 变换 法 
( differential transformation method , DTM ) 作为 一 种 半 
解析 方法 ,具有 求解 精度 和 计算 效率 高 的 优点 。 
文献 [20] 采 用 DTM 方法 对 微 板 振 动 的 控制 方程 进 
行 求解 并 验证 了 结果 的 准确 性 。 文 献 121 |] 考虑 梁 
表面 影响 ,利用 DTM 方法 求解 了 温度 场 和 电磁 场 作 
用 下 的 弹性 约束 压 电 纳米 梁 的 固有 频率 ,并 验证 了 
DTM 方法 相对 于 其 他 经 典 算法 具有 更 快 的 收敛 速 
度 。 文 献 [22 ] 采 用 DTM 方法 求解 流 固 耦合 问题 ,并 
将 所 得 结果 与 有 限 单元 法 进行 比较 ,证 明了 DTM 方 
法 的 计算 效率 更 高 。 因 此 ,DTM 是 一 种 有 效 的 分 析 
流体 输送 管道 振动 的 方法 。 对 于 流 固 看 合 振动 问 
题 , 文 献 多 利用 DTM 方法 研究 不 可 压缩 无 旋 无 黏 
流体 [622231。 

本 研究 考虑 流体 的 黏 弹性 ,应 用 DTM 方法 来 研 
究 输 流 微 管 的 振 劲 特性 ,通过 数值 求解 分 析 流 体 黏 
性 、 微 管材 料 的 内 豪 特 征 尺 寸 参数 和 预 应 力 对 前 三 
阶 固有 频率 的 影响 。 与 以 往 的 输 流 微 纳米 管 的 动力 
学 研究 相 比 ,有 如 下 创新 点 :首先 ,以 往 文献 多 采用 
非 局 部 理论 研究 微 纳米 管 特别 是 输 流 纳 米 管 的 动力 
学 特性 ,本 研究 利用 修正 偶 应 力 理论 结合 微分 变换 
法 研究 输 流 微 管 的 振动 及 其 临界 流速 ;其 次 ,考虑 实 
际 流体 黏 性 时 ,流体 在 微 管内 状态 为 层 流 ,通过 计算 
最 大 层 流 速度 ,明确 本 研究 流体 速度 范围 , 相 比 于 义 
速 假 设 , 本 研究 的 动力 响应 结果 更 加 接近 工程 实际 。 
此 外 ,考虑 到 输 流 管 一 般 会 承受 轴 向 力作 用 ,本 研究 
探讨 微 管 轴 向 预 应 力 对 流体 临界 流速 的 影响 ,因为 
临界 流速 与 输 流 微 管 的 输送 效率 密切 相关 。 


1 物理 模型 


建立 如 图 1 所 示 输 流 微 管道 坐标 系 ,不 可 压缩 
流体 速度 v 沿 x 轴 正 向 , 梁 在 两 支撑 间 的 长 度 为 工 ， 
宽 为 5, 厚 为 h。 微 管 x 方向 预 应 力 为 co, 弹 性 模 量 
为 ,密度 为 p, ,流体 密度 为 pi, 微 管内 径 为 Ri ,外 径 
为 R, ,流体 截面 积 为 4 , 微 管 截面 积 为 A0 

根据 YANG 等 1 提出 的 修正 偶 应 力 理 论 , 细 长 
微 管 采用 Euler 梁 模型 , 则 微 管 的 应 变 能 变 分 为 
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其 中 :! AA EH ERER SE SG w 为 沿 z 轴 的 横 
向 位 移 ;C = E/2(1 +v) 为 拉 梅 常数 或 剪 切 横 量 ,这 
E v 为 泊 松 比 。 式 (1) 考 虑 了 截面 内 的 内 力 在 微 段 
变形 上 所 做 的 虚 功 之 和 。 


图 1 输 流 微 管 


Fig.1 Schematic diagram of a fluid-conveying micro-tubule 


道 示意 图 


”在 定常 流动 中 ,流体 在 微 圆 管 中 层 流速 度 方 
jest 525) 


e £^ 


ApC R? -r) 
NET (2) 
X. Ri 为 微 圆 管 半径 ;r 为 径 向 坐标 ;为 流体 动力 


TR; Ap 为 管道 两 端 受到 的 流体 压力 差 。 
- ,由 流体 流量 9 = = f. 2mrvdr 以 及 流体 平均 速 


,可 以 得 到 微 圆 管 中 层 流速 度 分 布 方程 , 即 
um m (3) 
流体 质点 沿 z 轴 的 横向 加 速度 a HB 
LT 合成 定理 计算 , BU 
E dw. Fw. Fw, dw (4) 
dt ðt ðxðt ðx 


其 中 式 (4) 右边 分 别 代表 流体 沿 > ba 车 加 速度 P} 
氏 加 速度 和 向 心 加 速度 。 
Ca 


=- ff, (o. È m n E rw) dAde - 
[Lotti 
=- | f(a - KSu tp t 28v  dAd: - 


L, R 


JU. 


其 中 ,u SPIRI ERR 
上 式 考 虑 了 外 力 所 做 的 虚 功 之 和 ,包括 微 管 沿 
着 轴 向 和 横向 的 惯性 力 所 做 虚 功 ,也 包括 流体 与 微 


w w w, 
p2ar| T +w rum 十 zV Jdrõw)dx 
x 


(5) 
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管 相互 作用 惯性 力 、 离 心力 和 科 氏 力 的 影响 。 因 为 
本 研究 考虑 不 可 压缩 流体 ,因此 微 管 两 端的 压力 差 
不 做 功 。 同 时 , 微 管 受到 的 预 应 力 远 大 于 黏 性 力 , 因 


此 黏 性 力 虚 功 可 以 忽略 。 
将 式 (1) 和 式 (5) 代 入 虚 功 原理 表达 式 
ôU = êV (6) 
经 分 部 积分 后 得 到 


g ; ; 
EGAN 4 
^ Ox 


L ow - diu -, Pw 
Í (oa z + Mv Nv 2E Owdx + 
0 ðt oxot OX” 
Fw aw Fw 
A, -一 一 一 一 -04 — 
| vt ac P» ax? at 0 9/15 acd J 
gogo -— 
P x^ Ox 9 
CE7+CL4)| ， *pl —ów| + 
Ó w 0 0 
一 0w 
Ox 6 
ðw 4 E 
ooA—6w -0 (7) 
x 0 
其 中 


RI 4 
M = | pst ea Ri - p)dr, 


N= NL. (scr B E 
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由 式 (7) 可 得 到 控制 方程 , 即 
(EI + GPA,) A: 十 


Fw 
; + Mv v A 
ðt Oxot Ox 


-a2 -a2 
(oa o e Nos 


uu co ee (8) 
"^a? "7 ax ox 
两 端 简 支 边界 条 件 为 , 即 


w| -0 ^7 w 


= 0， 


Fw 
pa 
设 控制 方程 式 (8 ) 的 解 为 
w(x,t) = WGoe? (10) 
其 中 :W(x) 是 模 态 函数 ;w 是 固有 和 角 频 率 。 将 式 
(10) 代 入 式 (8) 得 到 


(pA W + psAsW - 


_ ðw 


=0 (9) 


x=0 ax? 


x=L 


2 i 
1. Toa? bm Yua 
Ox" Óx 


42 


W 
CEI + GPA, uy c0 (11) 


对 式 (11) 进 


PE 
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L ax Ox 
(EI + Gl A,) ow n (N? - TA) Fw -0 
p Ox L Ox 
(12) 
其 中 w = T x 一 rA 


2 DTM 应 用 


DTM 是 一 种 基于 泰勒 级 数 展开 的 用 于 求解 复 
杂 微 分 方程 的 变换 方法 , 它 使 用 充分 可 微 的 多 项 式 
形式 作为 精确 解 的 近似 ,使 计算 的 效率 大 大 提高 。 
根据 文献 127 ] ,微分 变换 和 微分 逆 变 换 分 别 表示 为 


inue S (EX (13) 
ie X (x - x)! YLk] (14) 


OK1 列 出 了 在 文章 分 析 中 常用 的 DTM 基本 
agm, 
R1 DTM 常用 基本 运算 
Tab.1 Basic mathematics operations of DTM 
wx) = y(x) £2(%) 


转化 后 的 函数 
W(k) = Y(k) + Z(k) 


(s) = Aya) WO) = AY(A) 
f x) = LL W(k) = (k+1)Y(k +1) 
gen aimo, W(k) = (k*1)(k *2)--(k * n)Y(Ck 9 n) 


9^ 


i L 
w(x) = y(x)z(x) W(k) = MY(ODZ(k - D) 
{=ü 


l,k=m 
w(x) = x" W(k) 2ó6(k-m) = : 
(x) (+) ( ) SENE 

对 式 (12) 微 分 变换 ,得 
WC +4) = 
F (pA LW + pA,LWCD - 
2 
+ 


PL -EDGO -2»)WCE +2) 


Mock + DWG «Do + 


OV! - oA) 
L——r DG «2WG €2) 


- CEI + GÊA,) 
~p & -rDO +2)(k +3)(k +4) 


(15) 
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类 似 地 ,边界 条 件 式 (9) 可 转化 为 


W(0) = W(2) =0 (16) 
Y wa) =o (17) 
Y kg DW) 20 (18) 


将 式 (15) 代 入 式 (17) 和 (18) ,所 得 两 个 方程 仅 
fur WCp I Wo 两 个 未 知 量 , 因 此 可 以 整理 成 如 
FABPAJE XX , HII 

Yu Ya] [W(1) 

网 ror nd 
其 中 , yy 与 微 管 系统 参数 和 固有 频率 o 相关 ,为 使 
式 (19) 有 非 零 解 ,需要 
Yu Yn 


=0 (20) 


Yu Yn 
由 式 (20) 可 以 解 得 各 阶 固 有 频率 ,图 2 给 出 了 
本 研究 计算 过 程 的 流程 图 。 


Wi70,w,, 76, 0,Wis=C, 
5,7€,5,76c, k-0,0,^ ?—0,6-0.001 


Wa Wk+4) 


mM=(k+4)(k+3)W, 
n=(k+4)(k+3Nk+2)W,., 


Sms tm, sy-s;tn 


Extract the coefficient of c, and c, in the s, expression, y,, and y,, 
Extract the coefficient of c, and c, in the s, expression, y,, and y, 


output w” 


图 2 微分 变换 法 计算 流程 简 图 


Fig.2 Calculation flow chart of differential 


transformation method 
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3 数值 算 例 


采用 环 氧 树脂 材料 ,弹性 模 量 E = 1. 44 GPa, 1H 
PSHE v 20.38,1 217.6 uum" 。 微 管内 径 尺 29.65 um, 
外 径 忆 =10 m. Ahi Hd SU VU 
小 的 情况 下 ,雷诺 数 计算 公式 为 
_ 2Ripw 
EE" 

对 微 管 Re «2.320 时 为 层 流 ,由 此 得 到 流体 在 微 
管 中 层 流 的 最 大 平均 速度 , 即 


- 2 320n 
v < ( 22 ) 
2R,p, 


Re (21) 
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WEWERE u =0.1 Pa - s, 则 由 式 (22) 
计算 最 大 层 流 平均 速度 为 15 026 m/s, 

表 2 将 DTM 方法 所 得 结果 与 有 限 差分 法 (finite 
differential method , FDM ) 所 得 结果 进行 了 对 比 ,其 中 
FDM 采用 100 个 计算 节点 ,两 种 方法 下 结果 误差 不 
超过 0.08% 。 图 3 给 出 了 前 三 阶 固 有 频率 与 流体 平 
均 速 度 的 关系 ,并 给 出 了 无 黏 性 流体 在 相同 平均 流 
速 条 件 下 的 固有 频率 。 输 流 微 管 的 振动 表现 出 流速 
越 大 固有 频率 越 小 的 特性 。 在 微 管内 流体 平均 速度 
相同 即 流量 相同 的 条 件 下 ,流体 黏 性 使 得 微 管 前 三 
阶 固有 频率 偏 低 , 且 平 均 速度 越 大 两 者 之 间 的 偏差 
越 大 。 


表 2 DTM 5 FDM 结果 对 比 
Tab.2 Natural frequency of tubules with different methods 


FDM DTM 
Mode 1/ Mode 2/ Mode 3/ Mode 1/ Mode 2/ Mode 3/ 
(108rad . s7!) (10frad - s7!) (105rad + s7!) (105rad - s7!) (10frad - s7!) (105rad - s7!) 


5 1.051 570 4.178 150 9.346 510 1.051 657 1.051 657 1.051 657 
25 0. 960 179 4. 107 850 9.282 110 0.960 263 0. 960 263 0. 960 263 
50 0.621 682 3.882 680 9.079 260 0.621 772 0.621 772 0.621 772 
hs NOT i Mode 1 
= Mode 2 
= Mode 3 
7:5 1 
个 | — tatik ` 
+ 一 一 itik 
S 
2 50 Ee 
- Mode 2 - 
-— € S 
E 2 
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图 3 微 管 固有 频率 与 流体 速度 的 关系 


Fig.3 Natural frequency of tubules vs. fluid velocity 


T 


值得 注意 的 是 ,流体 平均 速度 为 68 m/s 附近 
时 ,第 一 阶 固有 频率 为 零 ,这 一 流速 称 为 临界 速度 。 
此 时 输 流 微 管 失 稳 , 若 继续 增 大 流体 速度 则 流体 离 
心力 将 对 管道 结构 造成 破坏 ,这 一 失 稳 形式 通常 称 
为 动态 失 稳 ,在 微 管 系统 设计 时 需要 严格 避免 ;在 微 
管 两 端 固 支 的 条 件 下 , 微 管 系统 有 类 似 失 稳 形 式 ; 若 
微 管 悬臂 文 撑 , 根 据 现 有 理论 模型 可 以 断定 ,流速 增 
大 则 各 阶 固有 频率 降低 ,这 一 结论 与 宏观 尺度 管道 
相同 , 详 见 文献 [6]。 有 意思 的 是 ,悬臂 微 管 达到 临 


界 速度 后 ,继续 增 大 流速 不 会 造成 结构 破坏 ,因此 称 
为 静态 失 稳 。 限 于 篇 幅 , 固 支 和 悬臂 两 种 边界 条 件 
不 再 单独 给 出 计算 结果 。 

表 3 给 出 了 内 豪 特征 尺寸 对 前 三 阶 固有 频率 的 
影响 ,内 侵 特 征 尺 寸 越 大 则 前 三 阶 固 有 频率 越 大 。 
表 4 给 出 了 预 应 力 对 前 三 阶 固 有 频率 的 影响 ,在 预 
应 力 增 大 的 过 程 中 ,各 阶 频率 均 增 大 ,但 预 应 力 对 各 
阶 频率 的 影响 程度 不 同 , 对 第 一 阶 频率 影响 最 大 , 且 
预 应 力 影 响 随 阶 次 升 高 而 减 小 。 
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RI 微 管内 豪 特征 参数 与 固有 频率 的 关系 
Tab.3 Natural frequency of micro-tubules vs. the material 


length scale parameter 


Mode 1/ Mode 2/ Mode 3/ 
l/ um 
(105rad . s7!) (105rad . s7!) (10frad * s7!) 
10 0.369 886 2.775 460 6.551 480 
14 0.601 744 3.386 610 7.862 510 
17.6 0. 797 430 3.99] 850 9.176 930 
25 1.185 800 5.321 230 12.094 900 


表 4 ， 预 应 力 与 固有 频率 的 关系 


Tab.4 Natural frequency of tubules vs. prestress 


Mode 1/ Mode 2/ Mode 3/ 
g/ MPa 
(108rad . s7!) (10frad - s7!) (10frad - s7!) 
r 0.659 247 3.886 950 9.075 130 
0.785 852 3.982 430 9.167 730 
0. 895 067 4.075 670 9.259 390 
26.3596 4. 6396 1. 99% 


Nog - o, AD WCE +2) 
B 2 
— T +3)(k +4) 
> (23) 
1o, 代表 临界 流速 。 式 (23 ) 结合 边界 条 件 式 
D 和 (18) 可 以 解 得 各 阶 临界 流速 。 
.二 表 5 给 出 了 前 三 阶 临界 流速 与 预 应 力 的 关系 ， 
可 痪 预 应 力 增 大 使 得 各 阶 临界 流速 都 增 大 , 预 应 力 
对 第 一 阶 轩 有 频率 影响 最 大 , 预 应 力 的 影响 随 阶 次 
升 高 而 减 小 。 

R5 预 应 力 与 临界 流速 的 关系 


Tab.5 Critical speed versus prestress 


o/ MPa va/(m*s ^!) wg/(m*s^!)  vs/(m:s^!) 
-10 57.06 122.87 186. 60 
0 62. 84 125. 66 188.45 
10 68.12 128.38 190. 28 
增长 率 16. 24% 4.29% 1.93% 
4 结 it 


本 研究 基于 修正 偶 应 力 理论 和 Euler 梁 模 型 , 结 
合 虚 功 原理 推导 了 输 流 微 管 的 控制 方程 。 考 虑 流体 


投稿 网 站 :http://cjam. xjtu. edu. en 


JT I) 第 40 卷 

黏 性 ,给 出 了 流体 在 圆 形 微 管 中 的 速度 分 布 方程 。 
ANRE FERRA ENA AEM, E 
速度 分 布 不 均匀 。 在 平均 速度 相同 的 条 件 下 ,利用 
DTM 作为 求解 工具 ,探讨 了 速度 不 均 亦 即 流体 茜 性 
对 微 管 的 横向 振动 的 影响 ,并 将 数值 结果 与 有 限 差 
分 法 求解 结果 进行 了 对 比 。 讨 论 了 材料 内 豪 特 征 尺 
寸 和 预 应力 对 系统 振动 的 影响 。 求 解 了 前 三 阶 临界 
流速 ,分 析 了 预 应 力 对 临界 流速 的 影响 。 人 研究 表明 
黏 性 的 存在 使 得 前 三 阶 加 有 频率 趋 于 降低 ,流体 平 
yg JS BE ie d PESE WR S DRE EUBRAE S STER 
大 , 即 偶 应 力 效应 越 强 则 前 三 阶 轩 有 频率 越 大 。 随 
着 预 应 力 增 大 ,各 阶 固有 频率 和 各 阶 临界 流速 均 增 
大 , 阶 次 越 高 则 预 应 力 的 影响 越 小 。 
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